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ABSTRACT: Using Majorana Fermions to represent spins we re-examine the Kondo 
Lattice model for heavy fermions. The simplest decoupling procedure provides a realiza- 
tion of odd frequency superconductivity, with resonant pairing and surfaces of gap zeros. 
Spin and charge coherence factors vanish linearly with the energy on the Fermi surface, 
predicting a linear specific heat, but a T 3 NMR relaxation rate. Possible application to 
heavy fermions is suggested. 
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Though a decade and a half has passed since the discovery of the heavy fermion metals 
and superconductors 1,2 , many experimental anomalies remain. Whilst the basic theo- 
retical picture of resonantly scattered conduction electrons forming a highly renormalized 
f-band is not in question, certain experimental features fit awkwardly into the standard 
model . The underlying nature of the interactions , the nature of the pairing , and 
the excitation spectrum of heavy fermion insulators 10-13 are three areas of continuing 
uncertainty. 

Conventional approaches to heavy fermion physics represent the f-moments as fermions 
by enforcing a "Gutzwiller constraint" of unit occupancy rif = 1 at each site. This requires 
a projection of the physical Hilbert space of the local moments from the larger Hilbert space 
of pseudo-fermions: a task that is difficult to do exactly, and usually only treated on the 
average. In this letter we examine an approach to a simple Kondo lattice model for heavy 
fermions that avoids these difficulties. 

Various new features are predicted that differ qualitatively from the standard model 
of heavy fermion behavior; most notably a development of strong correlations between the 
spin and pair degrees of freedom, forming a ground state where the conduction electrons 
experience frequency dependent or "resonant" triplet pairing. The pairing fields actually 
diverge at low frequencies as the inverse frequency, providing a first stable realization of the 
phenomenon of odd frequency pairing originally considered by Berezinskii 14 ' 1 "'. For a wide 
range of conditions, including the presence of spin-orbit coupling, this theory predicts sur- 
faces of gapless excitations, and a linear specific heat that survives in the superconducting 
state. Unlike a conventional superconductor, the charge and spin coherence factors vanish 
on the pseudo-Fermi surface, giving rise to a T 3 NMR relaxation. In this scenario, the 
linear specific heat anomalies often observed in heavy fermion superconductors 16 might be 
interpreted as intrinsic. 

A key feature of our approach is the use of a special anticommuting representation of 
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spin one-half operators to describe the magnetic excitations within the low-lying crystal 
field doublets of the heavy fermion ions 17 . Recall that for individual S = 1/2 objects, 
the Pauli matrices are anticommuting variables {a a , 07,} = 2<5 a ^ and consequently can be 
treated as real or Majorana (at = a) Fermi fields. Their Fermi statistics alone guarantee 
that the spin operator S = — fa x a satisfies both the SU(2) algebra [S a , S b ] = ie a b c S c 
and the condition S 2 = 3/4. This feature can be generalized to many sites, introducing a 
set of three component anticommuting real vectors ffi at each site i, 

{vf,V b j } = S ij S ab (v] = V ] j a ) (a, b = 1,2, 3) (1) 

from which the spin operator at each site is constructed 

-* i 

S J = ~2^j x % ( 2 ) 
This "Majorana" representation of spin 1/2 operators has a long history 18 in particle 
physics. Loosely speaking, the Majorana fermions may be considered to be lattice gener- 
alizations of anticommuting Pauli operators ifj = ~^<?j- There is no constraint associated 
with this representation, for the spin algebra and the condition 5=1/2 are satisfied at 
each site, between all states of the Fock space 19 . In momentum space, the Bloch waves, 
fj^ = J2j rije~ lk ' R 3 behave as conventional complex fermions, but since 77^ = r/_-, the mo- 
mentum lies in one half of the Brillouin zone. Finally note that since there is no constraint, 
the trial ground-state energy obtained from a trial Hamiltonian is a strict variational upper 
bound on the true ground-state energy. 

Our basic model for a heavy fermion system is a spin 1/2 Kondo lattice model, with a 

— * 

single band interacting with local f-moments Sj in each unit cell. Our simplified Hamilto- 
nian is written 

H = H c + Y,H int [j] (3) 
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Here H c = describes the conduction band, and ^t- = (V't^, i s a conduc- 

tion electron spinor. The exchange interaction at each site j is written in a tight binding 
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representation as 

In a real heavy fermion system, we envisage that the indices would refer to the conserved 
pseudospin indices of the low lying magnetic manifold. We have suppressed both the 
momentum dependence and anisotropy of the coupling, using ia.(ff x if) = [if ■ a] 2 — | to 
simplify the interaction. 

We may now write the partition function as a path integral, Z = J p e ^( T ) rfr where 

^) = E%H+ E V^rffk + Hc + J^HinAj]- (4) 

k ke^BZ J 

Here we have factorized the interaction in terms of a fluctuating two-component spinor 

vij = (v*,v*) 

Hintij] = ■ Vj)Vj + V^0 ■ rfj)i>j + 2\Vj\ 2 /J (5) 

We are particularly interested in examining static mean field solutions where 

To gain insight into this mean field theory, let us integrate out the localized spin degrees of 
freedom, represented by the Majorana fermions. This introduces a resonant self-energy into 
the electron propagators, containing an isotropic component that builds the renormalized 
heavy fermion band and an anisotropic term, compactly represented by the effective action 

s c = j2 *uA^ +€ t +A(fj) h*> +Sa (7) 

{k,iiv n } 

where A (a;) = determines the strength of the resonant scattering. The anisotropic 
term S a is written in a tight binding basis as 

S a = - — + ■a}iJj, u) + ^,-uh?-<?-]^,w + c.c. } (8) 

{j, iu n } 



Here the triad of orthogonal unit vectors b = z^az, d = x + iy = z T [i<J2&] z define the 
orientation of the order parameter. The quantities 



B 3 M = ^M(y, = ^Md 3 (9) 



may be interpreted as resonant Weiss and triplet pairing fields, respectively. Unlike ear- 
lier realizations of odd frequency triplet pairing 14 ' 15 , here the odd frequency pairing field 
diverges at zero frequency, coupling spin and triplet pair degrees of freedom in one order 
parameter. 

To simplify further discussion, we consider the case of a bipartite lattice. Here, a stable 
mean field solution is obtained with a staggered order parameter, where for example b is 
constant, and d = e l Q' Rj d is staggered commensurately with Q = (tt, tt, tt). In this case 
the spinor Zj = e^'^/ 2 z , where z Q = ^o^- Writing the conduction electron spinors in 
terms of their four real components x^(/c) (A = 0, 1, 2, 3) 

V>j = ^{^j + i Xj -v}zo (10) 
the mean field Hamiltonian takes the simple form 

h mf = e {h4 ] 4 +iV b ] k^k- c - c ] +a k n k} ( n ) 



k&\BZ 



where 



k 



^M + ^M)-™] (12) 



is the number operator of the state k, written in the four component basis and 

6 k = 2 <y£ k+Q/2 ~ e -k+Q/2^ a k = 2 <y€ k+Q/2 + € -k+Q/2^ = ~^ ^ 

where the last equality holds only for a tight binding model. Let us begin by considering 
the special case of half filling, {ji = 0) , for in this case the Hamiltonian is diagonal in the 
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Majorana components A, with excitation energies 



%=| ± V(|) 2 + y2 « =1 ' 3 ) (14) 

corresponding to three hybridized gapful branches and a forth gapless Majorana mode 
formed from a component of the conduction band that does not mix with the local mo- 
ments. With one unpaired Majorana fermion per unit cell, the corresponding Fermi surface 
= spans precisely one half of the Brillouin zone: Vps/(2n)^ = \. This counting ar- 
gument guarantees that the gapless Fermi surface persists in the presence of particle hole 
asymmetry (fx 7^ 0) or a spin dependent kinetic energy associated with spin-orbit coupling. 

For our particular choice of z , the up electrons are "paired", whilst the "down" elec- 
trons are unpaired with a gapped excitation spectrum (Fig. 1.). In a Nambu notation, 
their propagators are 

_ ([(.-epAMfl+n)]- 1 (a=T) 
G a (u,k) = { (15) 
{ [(u-e^-2A(u)}- 1 (a=i) 

where = — fxr^. The density of states for the "up" electrons is 

f?(l + "~ r) (M < t k) 

Pl(u>) = I 2V [F 4 /4 + ^ 2 ]2 ; , (16) 

[p (\<o\>T K ) 

v 2 

where Tr = j^fi^jjyjT^ is the indirect gap associated with the excitation spectrum (14). 

Unlike conventional pairing, the charge and spin coherence factors of these "Majo- 
rana" quasiparticles are strongly energy dependent. Near the Fermi surface, the gapless 
quasiparticle operators can be written as 

V = ^M^t + + ^ - z ^3k ( 17 ) 

-1 u 2 

where Z =l + -^-isa quasiparticle renormalization constant and the Bogoliubov 
coefficients are determined by v? + v 2 = 1, 



2 1 

U k=2 



u sgn(E'r) 

1+ k = \ (18) 



^ 2 + " 2 



Spin and charge coherence factors are then given by 

(k-\p^k+) ] z 
{k |cr-|/e ) J 

These quasiparticles thus form a pseudogap where spin and charge matrix elements vanish 
on the Fermi surface and grow linearly with energy, In the special particle-hole symmetric 
case (u. = 0), these coherence factors vanish throughout the gap, forming a neutral band 
of excitations that only conduct heat. Since the paramagnetic spin and charge response 
functions of the quasiparticle fluid are proportional to the square of these matrix elements, 
the corresponding local response functions to grow quadratically with energy 

%^) 2 (£) 2 <*> 

This unusual energy dependence of matrix elements permits this state to mimic one with 
constant coherence factors, but a linear density of states (line of gap zeros). 

We briefly list the main consequences of these results: 

i) A large quasiparticle thermal conductivity in the absence of a quasiparticle contribution 
to the thermopower and electrical conductivity. 

ii) Linear specific heat coefficient of magnitude 7 = g7n(l + A* 2 /^ 2 ) where 7 n is the 
linear specific heat in the absence of the local moments. As [i varies 7 can vary 
between values characteristic of a conventional metal, and values characteristic of a 
heavy fermion metal. 

iii) A T 3 component to the NMR relaxation rate superimposed upon an activated back- 
ground. 



Since the spin matrix elements (elS^e) = 0, the T 3 response is anisotropic and vanishes 
when the applied field is parallel to the b axis. 
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Finally, we should like to mention the collective properties of this state. Past studies 
of odd-frequency pairing have encountered a negative phase stiffness 20 . In our mean 
field theory, the phase has "coiled up" into a staggered configuration: this stabilizes the 
state and develops a positive phase stiffness. To compute the London response to a vector 
potential A, we replace — ► ^^_ e ^ T • The London Kernel Ad a fj = V 2 ^ ^FfA] is then 

e 2 v 2 T 



A = 



6 

K 



]TTr [g c ^K) 2 - S cT («)t 3 S cT («)t 3 ] (22) 



where the minus sign is a result of the staggered phase. Carrying out the energy integral, 
at T = 

Ne 2 [ Tk AM 2 
A = — — / duo =^ 3 (23) 

4 ™ Jo (A(u;) 2 + ^ 2 )§ 

2 

where we have set ^ = In the special case of fj, = 0, this integral simplifies to 

A = ( ) This small phase stiffness is consistent with the large coherence lengths 



= /U A of heavy fermion superconductors, and may be regarded as a consequence of a 
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condensation of "heavy electrons" with mass m* = m^. 

Macroscopic properties of the paired state are governed by slow rotations of the or- 
der parameter. In the absence of anisotropy, the long-wavelength action is fully spin- 
rotationally invariant, given by a Uncharge x SU(2) sp i n sigma model. Despite the phase 
stiffness, it can not support a super current without anisotropy, for the vacuum is not topo- 
logically against the creation of vortices: uniform gradients of the phase can continuously 
develop to cancel an externally applied vector potential by twisting the order parame- 
ter "into the third dimension" 21 . Anisotropy aligns the order parameter with the crystal 
axes, lowering the symmetry to a ^ r (l) c / iarffe X — Y model, where a supercurrent is topo- 
logically stable. In the special case of half filling the order parameter can be uniformly 
rotated in particle-hole space ^ J* ^ — > g ^ J* ^ where g is an SU(2) matrix. Now the 
long-wavelength action possess an SU(2) c i iarge symmetry, and again, does not support a 
persistent current. Since the half filled state has a gap to both spin and charge excitations, 
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this suggests this state will be insulating: a "superconductor" with a vanishing critical 
current and upper-critical field. 

Several questions of a technical nature surround our simple mean field theory. One 
performance benchmark of our mft is provided by the one impurity Kondo model: here we 
may compare its performance with the exact results and the well established large iV mean 
field theory (mft) 3 . For this model, the Majorana mft correctly yields a unitary ^ phase 
shift for the scattered electrons; it also predicts an enhanced isotropic susceptibility and 
linear specific heat: features consistent with the Fermi liquid fixed point. The mean field 
Wilson ratio * = 8/3 compares more favorably with the exact value ^ = 2 than the large 
N mft, where ( ^ ) = 1. As in the large N approach, RPA fluctuations in the phase 

variables develop power-law correlations the mean field order parameter, but here there is 
no local continuous gauge symmetry so power-law correlations are physically manifested 
as long-time correlations of the spin-charge operator M(t) = S(t)p(t). Verification of 
such correlations in the one impurity model would provide an independent test of our 
technique. Beyond the one impurity model, it remains to be seen whether our approach 
can also recover the normal phase by a careful treatment of these fluctuations. 

Experimentally, the strong frequency dependence of coherence factors in our theory 
may help reconcile the observation of large linear specific heats and thermal conductivities 
in heavy fermion superconductors with the consistent absence of a corresponding Kor- 
ringa NMR relaxation normally associated with gap less superconductivity. There is also 
a possible link with Kondo insulators u > 12 . In CeNiSn, a T 3 NMR relaxation rate and 
pseudogap have been observed 12 , features consistent with incipient odd pairing. It would 
be interesting to measure and compare the thermal conductivity of CeNiSn with LaNiSn, 
to check for anisotropies in the T 3 NMR response and test for a possible proximity effect 
with other heavy fermion superconductors. 

In conclusion, we have examined an alternative treatment of the Kondo lattice model 
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for heavy fermions that uses a Majorana representation of the spins. Our theory predicts 
a low temperature ground-state with odd frequency triplet pairing and surfaces of gapless 
neutral excitations. Spin and charge coherence factors vanish on the Fermi surface, pre- 
dicting an intrinsic thermal conductivity and linear specific heat that coexist-exist with a 
T NMR relaxation rate. Independently of these issues, it provides a first stable realization 
of Berezinskii's odd frequency pairing. 
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FIGURE CAPTIONS 



Fig. 1. Excitation spectrum of mean field theory for \i = 0, showing three- fold degenerate 
gapped excitations and a gapless Majorana band. Inset, conduction electron density of 
states for up (f) and down (|) electrons. 
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